
LOAN COPY: RETURN TO 
AFWL (WLIL-2) 
K1RTLAND AFB, N MEX 


ANALYSIS OF THREE-FLUID, 

CROSSFLOW HEAT EXCHANGERS 

A J.yf'Tfe. 

by Noel C Willis, Jr. 

Manned Spacecraft Center f f/p 

Houston, Texas * " v " ^ 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


MAY 1968 


TECH LIBRARY KAFB, NM 



tech library kafb, nm 



□ 01,8425 


ANALYSIS OF THREE-FLUID, CROSSFLOW HEAT EXCHANGERS 


By Noel C. Willis, Jr. 

Manned Spacecraft Center 
Houston, Texas 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



ABSTRACT 


The detailed behavior of three -fluid, crossflow 
heat exchangers has been investigated. The equations 
governing the two-dimensional temperature distribu- 
tions of the three fluids have been derived and nondi- 
mentionalized. Performance characteristics have 
been determined for a wide range of operating param- 
eters for single -pass heat exchangers. The perform- 
ance of two -pass heat exchangers for both cocurrent 
and countercurrent flow has been studied for selected 
operating conditions. Results have been presented 
graphically in terms of the temperature effectiveness 
of the two outer fluids as functions of heat -exchanger 
size for sets of fixed operating conditions. Nondi- 
mensional operating parameters have been defined 
which allow an efficient presentation of the large vol- 
ume of performance data required to representa prac- 
tical range of operating conditions. Sample problems 
are included to illustrate the use of the performance 
graphs for design applications. 
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ANALYSIS OF THREE- FLUID, CROSSFLOW HEAT EXCHANGERS 


By Noel C. Willis, Jr. 
Manned Spacecraft Center 


SUMMARY 


The detailed behavior of three-fluid, crossflow heat exchangers has been investi- 
gated. The equations governing two-dimensional temperature distributions of the three 
fluids have been derived and nondimensionalized. The performance characteristics 
have been determined for a wide range of operating parameters for single-pass heat 
exchangers. Performance of two-pass heat exchangers for both cocurrent and counter- 
current flow has been studied for selected operating conditions. Results have been pre- 
sented graphically in terms of the temperature effectiveness of the two outer fluids as 
functions of heat-exchanger size for sets of fixed operating conditions. Interpolation 
techniques have been used to obtain performance data for intermediate values. Nondi- 
mensionalized operating parameters have been defined which allow an efficient presen- 
tation of the large volume of performance data required to represent a practical range 
of operating conditions. Sample problems are included to illustrate the use of the per- 
formance curves and the interpolation techniques. 

An expression for overall effectiveness has been derived which compares the heat 
transferred by a particular exchanger with that transferred by one of infinite size. Iso- 
lated cases, corresponding to poor design, are cited in which the overall effectiveness 
may be greater than unity. This effect emphasizes the importance of using the temper- 
ature effectiveness of the two outer fluids as the primary design variables and the over- 
all effectiveness as an auxiliary parameter. 

A computer program has been developed for the study of both single- and 
multiple-pass heat exchangers. Output options are available for detailed studies of 
temperature distributions within a particular exchanger and for generation of perform- 
ance data for a large number of heat exchangers. 


INTRODUCTION 


Considerable effort has been expended in previous investigations to define per- 
formance characteristics of heat exchangers involving energy transfer between two 
fluids. Now that industrial processes have been developed which require simultaneous 
heat exchange between more than two fluids, analytical techniques are needed to de- 
scribe the performance of multifluid heat exchangers. One example of such a process 
is the large-scale production of oxygen in an air separation plant which requires heat 
exchange between oxygen, nitrogen, and air at very low temperatures. There are also 



possibilities for combining several separate two-fluid heat-exchanging operations more 
economically in a single multifluid arrangement. 

Several investigators (refs. 1 to 5) have pursued the problem of multifluid heat 
exchangers in parallel or counterflow in which only one physical dimension of the ex- 
changer is considered. The purpose of this investigation is the detailed study of three- 
fluid heat exchangers in crossflow as a completely two-dimensional problem. 

In this study, the performance of three-fluid crossflow heat exchangers is deter- 
mined and presented graphically in terms of the temperature effectiveness of two of the 
fluids referred to the third fluid. The effectiveness is determined as a function of heat- 
exchanger size for sets of fixed operating conditions. The introduction of nondimen- 
sional operating variables reduces the volume of data required to represent a practical 
range of operating conditions. The number of boundary conditions for the temperatures 
is reduced from three to one by the introduction of a nondimensional inlet temperature 
parameter. 

An expression for overall effectiveness is derived which compares the perform- 
ance of a heat exchanger to that of an infinitely large exchanger operating at the same 
conditions. A study of the two-dimensional temperature distributions reveals circum- 
stances in which the overall effectiveness may be greater than unity. This result im- 
plies that in a three-fluid, crossflow heat exchanger, total heat transfer is not always 
maximized by increasing the size of the exchanger. 

Effectiveness factors are determined for a wide range of operating parameters 
for single-pass, three-fluid heat exchangers. Performance of multiple -pass, three- 
fluid heat exchangers for both cocurrent and countercurrent flow is studied for selected 
operating conditions. 

Sample problems are used to illustrate the application of the effectiveness curves 
to heat- exchanger design problems. Since some of the performance data can be ex- 
plained only in terms of the two-dimensional variation of the temperatures of each fluid, 
these problem solutions are also used to provide insight into the detailed behavior of 
the fluids within the heat exchangers. 

The basic differential equations for the spatial distribution of the temperatures of 
the three fluids were solved numerically using a digital computer. A program is avail- 
able for both single- and multiple-pass calculations. An automatic -integration step- 
size control was developed through the consideration of overall conservation of energy 
so that multiple cases may be run continuously with the optimum step size used for each 
individual case. Output options are available for a detailed study of spatial temperature 
distribution within the exchanger or for the determination of the overall performance 
using only the average exit temperatures and effectiveness values. 


SYMBOLS 


A 

u 1 0 x y 
1, 2 o J o 


"Vp, 1 

B 

u 2, 3 x o y o 

A 3 c p, 3 

p 

u l, 2 x o y o 


*2 c p, 2 

C P 

specific heat of fluid 

C p,j 

specific heat of fluid (j) 

D 

u 2, 3 x o y o 

A 2 c p, 2 

E 

Q a 

overall effectiveness, q- 

K 1 

m 1 c 1 

, capacity rate ratio for fluid (1) 

m 2 c p, 2 

K 3 

m„c o 

— — — , capacity rate ratio for fluid (3) 

m 2 C p, 2 

M 

largest value of the set (A, B, C, D) 

m 

mass flow rate of fluid 

NTUj 

number of transfer units of a heat exchanger referred to fluid (1) (equal 
to A) 

NTUg 

number of transfer units of a heat exchanger referred to fluid (3) (equal 
to B) 

Q 

Q 

nondimensional heat-transfer rate, . _ t n 

m 2 C p, 2(^1, i l 2,i) 

Q 

total heat-transfer rate in a heat exchanger of finite size 
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Q 0 


total heat-transfer rate in a counterflow heat exchanger of infinite size 


q j 


heat transferred to fluid (j) per unit time 


1 - t 2 i 

nondimensional temperature, t r — (all subscripts listed for t 

1, i 2, i 

are applied to T also) 
local temperature of fluid 


3, mix 


2, e, « 

t 0 (X 

2, max 


inlet mixing temperature of fluids (1) and (3), 


rmc . + m„c „t 0 . 

1 p, 1 1, i 3 p, 3 3, x 

rmc t + m„c _ 

1 p, 1 3 p, 3 


exit temperature of fluid (2) for an infinitely large heat exchanger 
=0) maximum value of t ^ along the y-axis for an infinitely large heat ex- 


t. 

3 

t. 

J> e 


t. 

j,em 
t. . 


changer, 


Ut l,i +t 3,i 


U + 1 

local temperature of fluid (j) 
exit temperature of fluid (j) 
average exit temperature of fluid (j) 
inlet temperature of fluid (j) 


U 


u 


conductance ratio, 


1,2 


u 


2,3 


u 


1,2 


overall conductance between fluids (1) and (2) 


u 


2, 3 


overall conductance between fluids (2) and (3) 


nondimensional coordinate of heat- exchanger surface, — 

o 


x 

x 


coordinate of heat- exchanger surface 
x- dimension of heat exchanger 
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Y nondimensional coordinate of heat-exchanger surface, — 

^o 

y coordinate of heat- exchanger surface 

y Q y-dimension of heat exchanger 


At. 

1 

0 1 

0 3 


inlet temperature parameter, 




temperature effectiveness of fluid (1), 


ti . - t- 
1, i 1, em 



temperature effectiveness of fluid (3), ^ em 

t 3, i ~ t 2, i 


Subscripts: 

1 fluid (1) 

2 fluid (2) 

3 fluid (3) 

c cold fluid 

e exit (temperature) 

em average exit (temperature) 

h hot fluid 

i' inlet 

ident identical order 

inv inverted order 

j fluid (j); general or typical reference to fluid (1), (2), or (3) 

m average value (of temperature normal to fluid flow direction (appendix A)) 

max maximum value 

min minimum value 
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mix 


inlet mixing (temperature) 

number of integration step in x direction (appendix A) 


N 


Superscripts: 
c corrected 

p predicted 


PROBLEM FORMULATION 


Derivation of the Governing Equations for Three -Fluid, 

Crossflow Heat Exchangers 

Figure 1 is a schematic representation of a single-pass, three-fluid, crossflow 
heat exchanger. Heat is transferred between the center fluid (2) and outer fluids (1) 
and (3); however, there is no heat directly transferred between the two outer fluids. 

The immediate objective is to determine the temperature distributions of the three 
fluids in the heat exchanger for a given size and operating condition. Once the temper- 
ature distributions are known, the heat transferred to each fluid may be calculated from 
average exit temperatures; subsequently, the performance of the heat exchanger may be 
evaluated. 

In this investigation, the following simplifying assumptions have been made to re- 
duce the complexity of the equations. 

1. Steady flow exists for the three fluids. 

2. Fluid properties are constant. This assumption is adequate where large dif- 
ferences do not exist between the temperatures of the fluid and the heat-transfer sur- 
face. In most cases, evaluation of properties at a mixed mean temperature is sufficient 
to correct for property variations. A discussion of the effect of temperature-dependent 
fluid properties may be found in reference 6. 

3. For a particular surface, the local conductance is constant and equal to the 
overall conductance. This assumption is consistent with steady flow and constant fluid 
properties if entrance effects are neglected. 

4. The heat exchanger is considered to be adiabatic. If the performance of any 
heat exchanger is degraded because of heat exchange with the surroundings, it can be 
insulated. Effects caused by such heat exchange are not of interest in this investigation. 

5. The effects of longitudinal and lateral conduction in the heat exchanger are 
neglected. These effects are important if large temperature gradients exist and will 
reduce heat-exchanger effectiveness. 
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6. There is no lateral mixing in any fluid. This behavior is closely approxi- 
mated in plate-fin heat exchangers and when flows are not baffled. This assumption 
preserves the two-dimensional character of the problem. If mixing were assumed, the 
analysis would be simplified considerably. 

Under the previous assumptions, the governing equations for the temperature dis- 
tributions in three-fluid, crossflow heat exchangers will be derived. Figure 2 repre- 
sents the heat-transfer surface between fluids (1) and (2). For a properly designed 
exchanger, outer fluids (1) and (3) will be either both hotter or both colder than the 
center fluid (2). For purposes of discussion during this derivation, the center fluid is 
arbitrarily assumed to be hotter than the two outer fluids; however, the resulting equa- 
tions are independent of the assumed temperature levels. 

In figure 2 the heat transferred per unit time into fluid (1) from fluid (2) across 
the elemental area dx dy is 


dq i = u i, 2^2 - t i) dx d y 


( 1 ) 


where t^ and t ^ are both functions of x and y. 

This expression may be equated to the energy increase per unit time of the ele- 
ment of fluid (1) between y and y + dy as it moves from x to x + dx, which is 



(2) 


X V 

Introducing nondimensional coordinates X = — and Y = — and equating the above ex- 

x o 

pressions for dq^, the resulting differential equation is 



u i, 2 x 0 y 0 

rhiC 
1 p, 1 


(*2 " * 1 ) 


(3) 


Following a similar procedure for fluid (3) gives 


^3 u 2, 3 x o y o 

3X ' *3 c p, 3 



(4) 
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A differential volume element of fluid (2) is bounded by two surfaces and is in 
thermal communication with both fluids (1) and (3). The energy transferred to this ele- 
ment of fluid (2) from the outer fluids is 


dq 2 = [ u l, 2 (‘l - ‘ 2 ) + U 2, 3(‘s - t 2)] dx dy 


( 5 ) 


As the elemental volume of fluid (2) between x and x + dx moves from y to y + dy, 
its thermal energy increase may be expressed as 



Equating these two expressions and introducing the nondimensional coordinates yields 


t„ u., ~x y u 0 „x y 

2 = 1, Tjyo , _ , _ 2, 3 o y o . _ 

Y m 2 c p, 2 U 2 ' m 2 c p, 2 1 3 2> 


Three simultaneous partial differential equations have been derived which define the 
temperatures of each of the three fluids as functions of both space coordinates x and 
y. The resulting equations are 



( 8 ) 
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where 


1 , 2V0 
*1 C P,1 


B 


U 2, 3 X o y o 
m 3 c p, 3 


C 


u i, 2V0 

m 2 C p, 2 


\ 


D = 


u_ 0 x y 
_ 2, 3 o J o 


m 2 c p, 2 


(9) 


and 


h=t. (x,y) (10) 


Under the assumptions of the problem, the nondimensional terms A, B, C, and D 
are constants which depend upon the heat- capacity rates of the fluids, the heat- 
exchanger dimensions, and the values of conductance, at the two heat-transfer surfaces. 
When A, B, C, and D are specified, along with the inlet temperatures of the three 
fluids, the equations may be solved for t^, tg, and t^ as functions of position 

throughout the heat exchanger. 

Reduction of the number of bound ary conditions . - To nondimensionalize the basic 
equations with respect to temperature t., all temperatures can be referred to the inlet 

temperature of the center fluid t 0 . by subtracting t 9 . from the temperature t. and 

dividing this quantity by the difference between the inlet temperatures of fluids (1) and 
(2), or by t 1 . - t 9 .. The resulting equations are 


3T 1 

"ax" = a ( t 2 " T i) 

3T 9 

ST = C ( T ! " T 2) + D ( T 3 " T 2) ( 


3T„ 

IX = B ( T 2 " T 3) 


( 11 ) 
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where 


T - !nV 
1 * 1,1 “* 2,1 


Tr> = 


*2 ~ *2, i 


2 t- . - t, . 

l,i 2,i 


T _ *3 " * 2,1 
3 * 1,1 “* 2,1 


( 12 ) 


The advantage of the above formulation is apparent when the boundary conditions are 
examined. They become 


Tj(X=0) = 1 
T 2 (Y=0) = 0 


t 3 (x=o) 




/ 


(13) 


The boundary conditions for any problem may be specified by a single quantity called 
the inlet temperature parameter 





(14) 


This parameter is the ratio of the temperature levels of the two outer fluids referred to 
the temperature level of the center fluid. The third boundary condition becomes 


3, i 


At. ' 
l 


The outer fluids may be numbered so that At varies between zero and 


unity. As an example, in the case of both outer fluids being hotter than the center fluid, 

fluid (1) is always the colder of the two hot fluids. When At. is unity, t. . equals 

1 h ^ 

tg .; and when At. is very small, t^ . is considerably greater than t^ .. 
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For example, consider two sets of inlet temperatures 


and 


t. . = 300° F ) 

t„ . = 100° F > 

2.1 

t„ . = 500° F 
o, 1 J 

t 1 . = 50° F ^ 

1.1 

*2,i = 0<> F 

t, . = 100” F 

Oj X J 


( 15 ) 


(16) 


In both cases, At equals 0. 5, and the problems are equivalent in the foregoing non- 
dimensional formulation. 

Now that the specification of the inlet temperatures has been reduced to a single 
parameter, any problem may be completely defined by the five quantities A, B, C, 

D, and At. . 

’ l 

Multiple-pass, three-fluid, crossflow heat exchangers . - The use of multiple 
passes has long been recognized as a possible method of improving heat- exchanger per- 
formance. The previous investigations reported in references 7 and 8 have considered 
the problem of multiple-pass, two-fluid, crossflow heat exchangers. 

While the previous derivation and discussion has been applied to determining the 
temperature distributions in a single-pass heat exchanger, the basic procedure can also 
be applied to each pass of a multiple-pass heat exchanger. The extension of the analy- 
sis to multiple passes involves only the proper specification of boundary conditions for 
each pass for the several possible flow arrangements. The configuration for a two- 
pass, three-fluid, crossflow heat exchanger is illustrated in figure 3. This particular 
arrangement is called countercurrent. A cocurrent arrangement is obtained by revers- 
ing the indicated direction of fluid (2). The flow detail for each pass is the same as 
that illustrated in figure 1. 
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The solution of the cocurrent case is straightforward; the outlet temperatures of 
the first pass become inlet temperatures of the second pass, and so on, for all subse- 
quent passes. The solution in any pass is independent of the solutions of all subsequent 
passes. The temperature distribution in any pass is found in the same manner as in the 
single-pass heat exchanger, with the exception that the initial temperature generally 
will not be constant along the inlet boundaries after the first pass. 

The problem is somewhat more difficult for the countercurrent heat exchanger. 
Figure 4 represents the mathematical configuration for this case. The difficulty arises 
because the initial temperature distributions are not completely known for either pass. 
Using the terminology in reference 7, there are "outer" boundaries where the initial 
conditions for the whole exchanger are given in "inner" boundaries which are effec- 
tively common to the two passes. The inlet temperature for fluid (2) in pass 1 along an 
inner boundary is dependent upon the solution in pass 2. Since it is evident that the 
solution in pass 2 is dependent in turn on pass 1, some iterative scheme is suggested 
using assumed distributions along inner boundaries. 

Following the basic scheme outlined in reference 7, an inlet value of is as- 
sumed for pass 1, and the resulting exit values of T^ and T^ are used as input for 
pass 2. The exit value of Tg for pass 2 is used in the second calculation for pass 1. 
This iterative procedure is continued until the average exit values of T^ for pass 2 (for 

consecutive iterations) differ by less than a given convergence criterion. For the cal- 
culations in this study, the convergence criterion selected was that the average values 
of T 2 for consecutive iterations would differ by less than 1 percent. Depending upon 

the particular case, the number of iterations required for convergence was from three 
to five. 

For more than two passes, values of would have to be assumed on all inner 

boundaries where needed, the number of such boundaries being one less than the number 
of passes. Numerical calculations in this investigation were confined to two-pass heat 
exchangers. 

Discussion of the num e ric al solution of the basic equations. - The equations have 
been solved numerically by a first-order, predictor-corrector integration scheme. A 
FORTRAN program has been developed for use with the IBM 7094 or Univac 1107/1108 
computers. The program can handle single-pass and multiple-pass calculations. An 
automatic step-size control, governed by the values of certain input parameters, was 
developed so that a large number of cases could be run continuously with the optimum 
step size used for each individual case. Output options are available which allow the 
user to study either detailed temperature distributions or overall performance charac- 
teristics. The details of the numerical procedure are described in appendix A, and the 
computer program is discussed in appendix B. 
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I 


RESULTS AND DISCUSSION 


The performance of a two-fluid heat exchanger may be expressed by a single de- 
pendent variable which is a function of two independent variables. For example, the 
overall effectiveness can be expressed as a function of the number of transfer units 
NTU of the heat exchanger and the capacity rate ratios of the two fluids. An investiga- 
tion of three-fluid heat exchangers requires consideration of two dependent variables 
and five independent variables. 


Nondimensional Independent Variables 

In the foregoing section entitled "Problem Formulation," it was noted that the five 
independent parameters which can be used to define a specific problem are A, B, C, 
D, and AL. However, these are not easily associated with the physical variables of 

the problem. 

To specify a particular problem in terms of quantities which are more useful to 
a designer, the quantities A, B, C, and D may be combined into a new set of param- 
eters which are more amenable to physical interpretation. The new parameters are 


u , o x y 

NTU 1 = A = - - 

1 m l c p, 1 


(17) 


TT _ C _ U 1,2 
~ D_U 2,3 


(18) 


„ m.,c . 
K = £ = 1 P> 1 

1 A m 2 c p,2 


(19) 


and 



m 3 C p, 3 

m 2 C p, 2 


( 20 ) 


The constant A is retained as the basic size parameter and is called NTU^, or 

the number of transfer units of the heat exchanger referred to fluid (1). This parameter 
represents the ability of the heat exchanger to change the temperature of fluid (1). A 
large value of NTU^ can result from a large physical size x Q y o , a high conductance 
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between fluids (1) and (2) u., 0 , and a small capacity rate for fluid (1) m.,c v All of 

these factors would make fluid (1) relatively easy to heat or cool. Therefore, the non- 
dimensional input parameter NTU ^ is a good representation of the size of the ex- 


changer. Since a similar parameter based on either fluid (2) or (3) also could have 
been defined, the choice of fluid (1) as a reference for size is arbitrary. For example, 

u 2 3 x o y o 

a size parameter based on fluid (3) could have been defined as NTUg = — — 


is equal to the quantity B. 


m„c „ 
3 p, 3 


which 


The parameters and Kg are nondimensional heat- capacity flow rates of the 

outer fluids (1) and (3) referred to the center fluid (2) and will be called capacity rate 
ratios. For a well-designed heat exchanger, the combined heat-capacity flow rates of 
the outer two fluids should not be significantly different from the capacity flow rate of 
the center fluid. This implies that + Kg should be near unity for proper design. 


The parameter U is called the conductance ratio, and indicates the relative abil- 
ity of fluids (1) and (3) to transfer heat to fluid (2). 

Problems may be specified now by the five independent parameters NTUj, K^, 
K 3 , U, and At.. 


Discussion of Nondimensional Dependent Variables 

The solution of the basic equations provides two-dimensional distributions of the 
temperatures of all three fluids throughout the heat exchangers. Some of the phenom- 
ena which occur in three-fluid heat exchangers can be explained only by a study of these 
detailed distributions. Particular examples will be discussed in later sections. How- 
ever, the designer is interested mainly in the overall performance characteristics of 
the heat exchanger, for example, the average exit temperatures of the fluids. 

The dependent variables chosen to represent the performance of three-fluid heat 
exchangers are the temperature effectivenesses of the two outer fluids. These varia- 
bles are defined by the following expressions. 


1 - - — ^ x 100 percent 

1 

> 

0o =f44 em j x 100 percent 

3 \ *3,1 “*2,1/ 


( 21 ) 
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The quantities t 1 and t„ are the average exit temperatures of fluids (1) and 
lj 6m Oj 6m 

(3). They are obtained by averaging the local values for the exit temperatures obtained 
from the two-dimensional numerical integration of the basic equations. 

The variables 0^ and 0g represent the degree to which the temperatures of the 

outer fluids have approached the inlet temperature of the center fluid when they leave 
the heat exchanger. The effectiveness 0^ or 0g will be 100 percent when the average 

exit temperature of fluid (1) or fluid (3) equals the inlet temperature of fluid (2). The 
effectiveness will be zero when there is no change in temperature. There are circum- 
stances when one of the temperature effectivenesses actually can be negative. It has 
been assumed that for proper design, the center fluid (2) will either heat or cool both 
outer fluids. Consider the case for which t 1 . and t„ . are both greater than t 9 . 

1} 1 Uj 1 1 

and for which the function of the heat exchanger is to cool fluids (1) and (3). If the inlet 
temperature and heat capacity rates of fluid (3) are considerably greater than those of 
fluid (1), heat will be transferred through fluid (2) to fluid (1) in a large part of the heat 
exchanger. Fluid (1) will leave the heat exchanger at a temperature above its inlet 
value, the opposite of the desired effect. Therefore, a negative temperature effective- 
ness indicates that a fluid intended to be cooled was actually heated, or vice versa. A 
designer must use 0^ and 0g to determine the effect of the heat exchanger on both 

fluids. 

An auxiliary dependent variable, the overall effectiveness, has been defined to 
compare the performance of a particular heat exchanger to one of infinite size. The 

Q a 

overall effectiveness is E = where Q is the total heat transferred by an ex- 

^oo a 

changer under fixed operating conditions and is the heat that would be transferred 

by an infinitely large counterflow heat exchanger operating at the same conditions. It 
will be shown in a succeeding section that the maximum heat transfer does not neces- 
sarily occur for an infinitely large heat exchanger and that with the previous definition, 
the overall effectiveness can actually be greater than unity. 


Method of Presentation of Results 

The independent variables 0^, 0g, and E are functions of the independent, di- 
mensionless exchanger parameters K^, Kg, NTU, U, and AF. Since much heat- 

exchanger design work involves sizing an exchanger for a particular application, the 
performance factors 0^, 0g, and E are presented as functions of NTU for fixed K^, 

Kg, and U with AF as a parameter. Results for single-pass exchangers are pre- 
sented in figures 5 to 31 and results for two-pass exchangers in figures 32 to 37. 


Single -Pass Results 

A range of values for the independent variables has been chosen to cover a real- 
istic spectrum of operating conditions for the single-pass calculations. The variation 
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of the inlet temperature parameter can be confined to the range 0 to 1 by appropriately 
numbering the fluids. For example, in the case for which t.. . and t„ . are greater 

X, 1 o, 1 

then t 0 ., At. will always be 1. 0 or less if the colder of the two outer fluids is desig- 

nated as fluid (1). The conductance ratio U may vary from 0 to °°; however, the selec- 
ted values of 0. 5, 1. 0, and 2. 0 should be sufficient to cover the range of practical 
interest. The heat-capacity rate ratios and Kg may also vary from 0 to how- 
ever, Kj + Kg must be reasonably close to unity for a balanced heat exchanger. 
Therefore, the values for K^ and Kg of 0. 25, 0. 5, and 1. 0 should adequately cover 
the range of interest. The variation of NTU^ from 0 to 7. 5 is also sufficient to cover 

the range of practical sizes. The multiple-run option of the computer program was 
used to determine performance factors for heat exchangers represented by all possible 
combinations of the following set of independent parameters. 


NTtL = 0. 25 
NTUj = 2.0 


NTUj = 0. 50 
NTUj = 3.0 
NTUj = 7. 5 
Kj = 0.25 
Kj = 0. 50 

Kj = 1. 00 

Kg = 0. 25 
Kg = 0. 50 

Kg = 1. 00 


NTUj = 1.0 N 

NTUj = 5. 0 > 

) 


\ 

> 

/ 

\ 

> 

/ 


( 22 ) 


(23) 


(24) 
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and 


U = 0. 50 N 
U = 1. 00 . > 

U = 2. 00 


( 25 ) 


At. = 0. 25 

i 

At. = 0. 75 

l 


At. = 0. 50 

l 

> (26) 

At. = 1. 00 


This set of parameters required 756 separate calculations of the two-dimensional var- 
iations of all three-fluid temperatures. The printout option which restricted the output 
to the overall performance factors 0^, 0g, and E was used. Automatic step-size 

control was used to obtain the required accuracy and to minimize computation time. 

In the previous set of independent parameters there are three values each for K^, 

Kg, and U; therefore, there are 27 resulting performance charts for single-pass heat 

exchangers. Even with this many charts, only discrete values of the independent pa- 
rameters are represented. The succeeding section entitled "Application of Perform- 
ance Curves for Design" presents an interpolation technique for investigating problems 
defined by intermediate values of these parameters. 

The general trends of the data presented in figures 5 to 31 are similar to most 
heat-exchanger performance data with some exceptions. In most cases, the effective- 
ness factors increase with size, rapidly at first, then more slowly tending to some 
upper limit. This is always true for E and 0 gj however, in some cases 0^ begins 

to decrease as size increases and even becomes negative at times (fig. 11, for exam- 
ple). This effect may be explained as follows. Assuming an original intent to cool 
fluids (1) and (3), fluid (1) will be the coolest of the hot fluids. When 0^ decreases as 

size increases, fluid (1) has been actually heated in some portion of the heat exchanger 
which has been added, thereby reducing the effectiveness of the exchanger with respect 
to fluid (1). This occurs whenever fluid (2), by virtue of receiving heat from fluid (3), 
has been heated locally to a temperature greater than that of fluid (1). This effect will 
be discussed in detail in the succeeding section entitled "Overall Effectiveness. " When- 
ever the overall effect of the heat exchanger has been the heating of fluid (1) and the 
original intent was to cool fluid (1), the temperature effectiveness 0^ will be negative. 

Other trends can be noted which may be interpreted in terms of physical varia- 
bles. The patterns of behavior of the performance factors are dependent upon the rela- 
tive size of the heat exchanger referred to fluids (1) and (3). When the size variable 
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NTU^ was previously defined as 


able NTUg could be defined as 


u i, 2V 0 

a i c p, i 

U 2, 3 X o y o 
m 3 C p, 3 


= A, it was also noted that a similar vari- 


B. These parameters are measures of the 


ability of the heat exchanger to heat or cool fluid (1) or (3) per degree of temperature 
difference between that of fluid ((1) or fluid (3)) and fluid (2). 


The key to interpreting the behavior of the performance factors is the relative 
magnitude of A and B. Different patterns are noted for A > B, A = B, and A < B. 

In terms of the variables used in the performance charts, these criteria become 

K 1 K i K 1 K 1 

< U, = U, and ^ - > U. From the definition of the previous parameters, ~~ - U 

K 3 K 3 K 3 o 


fm.c 


l^p, 1 m 2 c p, 2)_ u l, 2 


is equivalent to 


u 


2,3 _ U l, 2 


, m 0 c „ m„c „ 
3p,3 2 p, 2 


u r 


or 


sides by X Q y o gives 

K 1 

and for A < B, ■=?- > U. 

K 3 


2, 3 

u, „x y u„ „x y 
1, 2 o J o _ 2, 3 o J o 


m 3 C p, 3 


m,c 
1 p, 1 


Multiplying both 
K„ 


m,c , 

1 p, 1 


m„c „ 
3 p, 3 


or A = B; similarly, for A > B, 


K, 


<U, 


The performance curves which fall under the above classifications are 



K 1 

A > B (figs. 7, 9, 10 to 13, 19, 21, 22, and 31); = U, A = B (figs. 6, 8, 16, 18, 

3 

K 1 

20, 28, and 30); and ~ > u, A < B (figs. 5, 14, 15, 17, 23 to 27, and 29). In cases 

K 3 

where A = B, the ability of the exchanger to heat or cool fluids (1) and (3) depends only 
on the temperature differential with the center fluid. When the inlet temperatures t 1 . 

1, X 

and t„ . are equal (At. = 1), 0 1 and 0„ are identical. For smaller values of At., 

Of i i i o i 

the temperature differential between fluids (3) and (2) is greater than that between 
fluids (1) and (3). The result is that heat is more easily transferred to fluid (3), and 
is always greater than 0^. 


When A < B, heat is transferred more easily to fluid (3) than to fluid (1). There- 
fore, in all cases when A < B, 9^ will be higher than 0j, even for At. = 1. It should 

be recalled that fluid (3) has been numbered so that At. is never greater than unity. 


There are two effects to consider when A is greater than B, the relative size 
(NTU) and the temperature differential. The 'size effect tends to make 0^ greater than 

0 3 , When At. = 1, 0j will always be greater than 0^; however, as At decreases, 
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the temperature differential between fluid (3) and fluid (2) becomes great enough to 
overcome the effect of NTU^ being greater than NTUg. The net result is that as ac- 
tual size increases, 9 g eventually becomes greater than 9 ^, when AU is less than 
1 . 0 . 


Multiple- Pass Results 

Numerical results have been obtained for the performance of representative con- 
figurations of two-pass, three-fluid, crossflow heat exchangers. 

For both the cocurrent and countercurrent cases there are several possibilities 
for the behavior of the fluids in the elbow sections of multiple-pass exchangers. A 
fluid may be completely mixed so that it enters one pass at a constant temperature — 
the average of the exit distribution from the other pass. The fluid may be completely 
unmixed in the elbow and approach the next pass with the identical temperature distri- 
bution with which it left the previous pass. Another possible condition would be no 
mixing in the elbow with a flow arrangement to invert the fluid prior to entering the next 
pass. Thus, for each fluid the following possibilities for the behavior in the elbow can 
be considered. 


1 . 

Mixed 


2. 

Unmixed, 

identical order 

3. 

Unmixed, 

inverted order 


All possible combinations would give 27 different cases to consider for each set 
of input variables K^, Kg, and U. To restrict the amount of data presented and still 

obtain an insight into the performance of three-fluid, multiple-pass, crossflow heat ex- 
changers, numerical results have been obtained for a countercurrent exchanger under 
the conditions = 0. 5, Kg = 0. 5, and U = 0. 5, 1. 0, and 2. 0 with all fluids mixed in 

the elbows. To provide comparisons for some of the other possibilities, the following 
cases have been considered for K n = 0. 5, K„ = 0. 5, and U = 1. 0: 


1. Cocurrent mixed 

2. Countercurrent, unmixed identical 

3. Countercurrent, unmixed inverted 

Figures 32 to 37 are the performance curves for the preceding cases. 

The size parameter NTU^ for multiple-pass cases is NA where N is the num- 
ber of passes and A, as previously defined, is ^ ^ ° ° . All other parameters are 

m lp, 1 

defined in the same manner as in the single-pass analysis. 
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To compare the performance of two-pass, countercurrent heat exchangers with 
single-pass exchangers operating at the same conditions, one can compare figures 17, 
18, and 19 with figures 32, 33, and 34, respectively. Although all of the effectiveness 
terms are increased in the two-pass, countercurrent case, the most significant in- 
creases occur in 0^ for the smaller values of Ah, particularly 0. 25. To obtain a 

direct comparison, the results for K^ = 0. 5, Kg = 0. 5, U = 1. 0, and At. = 0. 25 are 

presented in figure 38. 

It was noted previously in the case of the single-pass heat exchanger that for 
small At., the possibility existed for fluid (1) to be heated in some parts of the ex- 
changer, although the original intent was to cool fluids (1) and (3). This effect was the 
result of fluid (2) being heated by fluid (3) to a level that exceeded the local temperature 
of fluid (1). This tendency is decreased by the use of multiple passes in a counter- 
current arrangement. Fluid (1) may still be heated in pass 1; but in pass 2, fluid (3) 
has been cooled sufficiently so that it no longer heats fluid (2) above the level of 
fluid (1). The result is that fluid (1) is well cooled in pass 2; therefore, a significant 
increase in 0. is obtained for small At.. 

1 l 

Figure 35 presents the performance factors for a two-pass cocurrent arrange- 
ment for Kj = 0. 50, Kg = 0. 50, and U = 1. 0, with mixed flow in the elbows. This 

configuration is considerably less effective than the single-pass exchanger for similar 
conditions (fig. 18). The effect of the second pass was to reheat fluid (1) after it had 
been cooled in pass 1. The effectiveness decreased with an NTUj greater than 2. 0. 

These results for a specific case should not categorically condemn the multiple-pass, 
cocurrent arrangement, but they definitely illustrate the potential problem of reheating 
(or recooling) associated with this configuration. 

Figure 36 represents the performance factors for a two-pass, countercurrent ar- 
rangement for Kj = 0. 50, Kg = 0. 50, and U = 1. 0 for unmixed, identical- order flow 

in the elbows. Figure 37 represents a similar case for inverted flow; and figure 33, 
for mixed flow. Inspection of the curves indicates that all performance factors are 
highest for inverted order and lowest for identical, with mixed being slightly above 
identical. The differences are most pronounced for 0^, when Ah = 0. 25. Differences 

in overall effectiveness are slight; for example, for NTUj = 4. 0, E inv = 78 percent, 

E m ix = 76 percent, and = 75 percent. However, for Ah = 0. 25, 

0, . =45 percent, 0, . = 39 percent, and 9 . ., , = 37 percent for NTU, = 4.0. 

l,xnv 1, mix ^ ’ l,ident 1 

Again, general conclusions may not be drawn from this specific case. The case does 
indicate that the differences between the three possibilities for flow in the elbow are 
worth considering in design and may significantly affect some of the performance fac- 
tors of the exchanger. 

Figure 39 is an illustration of the convergence of T 9 . to pass 1 for the case of 

a, 1 

countercurrent flow with the inverted order Kj = 0. 5, Kg = 0. 5, U = 1 . 0, 

NTUj = 3. 0, and Ah = 0. 25. Four iterations were required. The initial estimate for 
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T 0 . to pass 1 was automated for the computation of multiple cases. The initial value 
1 

was assumed to be one- half of the mixing temperature of fluids (1) and (3) for each 
case. 


OVERALL EFFECTIVENESS OF THREE-FLUID, 
CROSSFLOW HEAT EXCHANGERS 


It is useful to define a parameter which compares the performance of a particular 
heat exchanger to a counterflow heat exchanger of infinite heat-transfer area operating 
at the same conditions. While the temperature effectivenesses 8^ and 0^ are of pri- 
mary interest to the designer, the overall effectiveness provides additional insight into 
heat- exchanger performance. The following discussion presents some interesting 
properties of three-fluid, crossflow heat exchangers that are most easily recognized 
and explained in terms of overall effectiveness. The overall effectiveness has been 

Q 

previously defined as E = where Q is the heat transferred by a particular ex- 

a 

changer and is the heat transferred by an infinitely large counterflow exchanger 
operating at the same conditions. The heat transfer Q is obtained from the solution 

* cl 

of the basic equations and may be expressed as either 


Q a " A 2 c p, sfa, 


em 


- t 


2,i) 


(27) 


or 


Q = m-c 1 ft 1 . - t 1 ) + m„c „ft„ . - t„ ) 
a 1 p, 1\1, i 1, emj 3 p, 3\3, i 3, em/ 


(28) 


In nondimensional form the two previous expressions become 




Q„ 


a m 


2 C p, 2( t l,i " l 2,i) ^ 


em 


(29) 


or 


Q a =K 


- (l - T, ) + K„ - T„ 

1( 1, emj 3 lAti 3,em 


(30) 
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Since there is some numerical inaccuracy in the computer program, Q is equated to 

CL 

the average of these two quantities. 

In deriving an expression for for a three-fluid heat exchanger, it is instruc- 
tive to consider a similar problem for a two-fluid heat exchanger, as illustrated in fig- 
ure 40. In the two-fluid, counterflow case, the exit temperature of the fluid with the 
smaller capacity rate will approach the inlet temperature of the fluid with the larger 
capacity rate as the heat-transfer area becomes infinite. If the hot fluid is denoted by 
the subscript h and the cold fluid by the subscript c, then for m^c^ h < m^c^ c , 

t, = t and for m,c , > m c , t = t, .. Therefore, the heat transfer for 
Ti, e c, x 5 h p, h c p, c’ c, e n, x 

an exchanger of infinite area is 


Q = ( 

^■OO 


me 


V R 


>) . Oh,! 

'mxn ' ’ 



(31) 


Figure 41 depicts the analogous situation for three-fluid, counterflow heat ex- 
changers. In the case for which the capacity rate of the center fluid is greater than the 
sum of the capacity rates of the outer fluids, the exit temperatures of the outer fluids 
both approach the inlet temperature of the center fluid. 


For 


( 


m,c . + m„c 
1 p, 1 3 p, 3; 


< m 2 C p, 2’ 


ti ^ = t 0 _ = t„ . and 
1 , e, °° 3, e, “ 2 , 1 


Q= m.c ./t. . - t„ .W m„c „(t~ . - t 0 The equivalent nondimensional quan- 
1 P? 1 , 1 1 j o p, o\ o, 1 <s, 1 / 


tities for 


( K i + ■%) 


< 1 are 


and 


T l> e , 00 ~ T 3, e, °° ~ T 2, i _ 0 


Q. 


K r 




flrjC Jt, . - t„ .) - K 1 + At. 

2 p, 2 ( 1, 1 2 , 1 ; 1 


(32) 


(33) 


For the situation in which ^m^ 1 + rfl 3 C p 3 ) > rfl 2 C p 2’ *2 e a PP roaches a limiting 

value t« which lies somewhere between t. . and t„ .. This limiting value 

2, e , 00 1 , 1 3 , 1 

should correspond to a single, effective inlet temperature of fluids (1) and (3). This 

effective inlet temperature is assumed to be the mixing temperature of fluids ( 1 ) and 

(3) defined as 


t 2, e , 00 t l, 3, mix 


rn.c ,t. . + m,c „t„ . 
_ 1 p,l l,i 3 p,-3 3 , 1 

m.c . + m„c „ 

1 p, 1 3 p, 3 


(34) 
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so that 


Q. 


^ 2 %, 2( t 2,e,°° 


5,0 


(35)' 


The equivalent nondimensional expressions for + Kg) > 1 are 


and 


K r 


K, +■ 


2, e, ~ 


K l T l a i " K 3 T 3,i _ 1 At i 


K l +K 3 


K l +K 3 


(36) 


5 X 

k 3 + At 

W °° m 9 c 0 (t. . - t 9 .) l 2, e, °° K 1 

2 p, 2 v 1, i 2, i) 1 , i 

K„ + 1 




(37) 


K 1 

It should be noted that T„ depends only on the ratio and At. and is independent 

2 , 6 Kg 1 

of U. 


To understand the variations of the overall effectiveness expression for three- 
fluid, crossflow heat exchangers, it is necessary to discuss the individual fluid temper- 
ature distributions for certain limiting situations. For two-fluid heat exchangers, the 
effectiveness can never be greater than 100 percent; however, using the above defini- 
tion, there are isolated cases corresponding to poor design practice for which the over- 
all effectiveness of a three -fluid heat exchanger in crossflow can actually exceed 
100 percent. A specific sample will be used to explain the behavior of the fluids when 
E is greater than 100 percent and to show why this somewhat anomalous result corre- 
sponds to poor design practice. The case under consideration will be one of using the 
center fluid to cool the two outer fluids. 

K l + 

the configuration is such that T 9 is greater than T„ = T. „ . . Figure 42 

represents the performance of a three-fluid, crossflow heat exchanger in terms of 0^, 
0g, and E for = 2. 0, Kg = 0. 5, and U = 2.0. The value of E reaches a maxi- 
mum of 101. 5 percent at NTUj = 4 for At. = 0. 25 and begins to decrease toward 


Values of effectiveness greater than 100 percent occur when 


Kg) > 1 and 
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100 percent as NTU^ continues to increase. To understand why it is possible for 

T 0 to be greater than T 0 , consider the distribution of T g along X = 0 from 

Y = 0 to Y = 1. 0. Along this line, and Tg are constant at their inlet values 

T 1 . and T„ . . Along the narrow strip at X = 0, fluids (1) and (3) act as infinite 

1, 1 o,i 

sources between which fluid (2) must flow. If the heat exchanger becomes infinitely 
large, there is a maximum temperature which fluid (2) may reach along the y-axis. 
Consider the case for T„ . > T.. . > T 0 .. 

1 lj 1 Ci} 1 

The maximum value of T 0 along the y-axis T„ (X=0) is reached when the 

ci Li } max 

heat flow rate into fluid (2) from fluid (3) is equal to heat flow rate from fluid (2) into 
fluid (1). This condition is expressed by 


u 


1, 2( T 2, max " T l, i) “ u 2, 3 ( T 3, i 


- T 


2, max) 


Solving for T 2 



max 


(X=0) 


UT 1 . + T 


U + 1 


3,i 


or since T- . = 1 and T„ . 

1, i 3, i 


1 

At 



max 


(X=0) 


U 


1 

At. 


l 

“ U + 1 


(38) 


(39) 


(40) 


The value T 0 is a function of At. and U, and is independent of K 1 and K„; 

ci j max i i o 

therefore, T„ is independent of T„ for fixed At.. If the value of T„ 

2, max r 2, e, 00 i 2, max 

is greater than Tg g /T^ g m ix) ’ effectiveness E may be greater than 
100 percent for sufficiently large values of NTU^ (the size parameter). 

Figure 43 illustrates the distribution of T^, Tg, and Tg along the y-axis for 
NTUj = 7. 5. Tj and Tg are constant, and Tg asymptotically approaches 2. For 
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(41) 


At. = 0. 25 and U = 2 
x 



max 


(X=0) 



U + 1 


2+4 

3 


= 2 


Figures 44 to 47 illustrate the temperature distributions as functions of Y for 
X = 0. 1, 0. 25, 0. 5 and 1. 0. As X becomes larger, all fluids approach a nondimen- 
sional temperature of 1. 6 at Y = 1. 0. This is more clearly illustrated in figure 48 in 
which all temperatures are plotted versus X at Y = 1. For these conditions, 


T 2, e, 00 ^1, 3, mix 


K 1 1 

K 3 +At i CT5 
K 2 

i r + 1 

K 3 


0. 25 


= 1.6 


0. 5 


+ 1 


(42) 


If the heat- exchanger size were made infinite, the exit temperature of fluid (2) 

would be T„ as defined previously, since the amount of fluid at a temperature 

z, e, 

greater than T„ would be negligible. For T 9 (X-0) greater than T„ , 

Uy Gj 00 max Zj 6j 

the overall effectiveness E will rise as NTU^ increases, will reach a maximum 
value, and will decrease asymptotically to 100 percent as NTUj becomes infinite. If 
Tg max (X=0) is * ess t* 13 - 11 T 2 e <»» the average exit temperature of fluid (2) approaches 
Tg e oo as a maximum, and E increases monotonically to 100 percent as NTU^ be- 
comes infinite. 

Examining the behavior of 0 0^, Tj, and Tg for this case, it can be seen 

that this behavior corresponds to poor design. The average temperature of fluid (1) is 
not affected and actually exceeds its inlet value at many points in the heat exchanger. 
Also, it can be seen that very little heat transfer is accomplished after X - 0. 5, indi- 
cating a highly oversized exchanger. 

This example should help emphasize the fact that while the overall effectiveness 
E is a very useful parameter, a designer should direct his attention to 6 ^ and 0g to 

analyze effectively the actual performance of the exchanger. 
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The following tables will assist the designer in quickly identifying the cases for 
which T 2j max (X=0) is greater than ^ ^ 


U + At 

T 8,m«»- 0 ) 'UTT 1 


2, max 


u 

At. = 0. 25 
1 

At. = 0. 50 
1 

At. = 0. 75 
1 

At. = 1.0 
1 

0. 5 

3.00 

1. 667 

1.22 

1.00 

1.0 

2. 50 

1. 50 

1. 167 

1.00 

2.0 

2. 00 

1. 33 

1. 11 

1.00 


_1 J_ 

K 3 + At. 

2, e , 00 K 

!T +1 

K 3 


K 1 


T 2,< 

00 

K 3 

At. = 0. 25 
1 

At. = 0. 50 
1 

At. = 0. 75 

X 

0.25 

3.40 

1.80 

1.267 

1.0 

2. 50 

1. 50 

1. 167 

2.0 

2. 00 

1.33 

1. Ill 

4.0 

1.60 

1.20 

1.067 


At. = 1. 0 
1 


Since T«, depends only on U and At., and T 0 depends only on the 

A) max x ci ) Gj 00 

K i 

atio rr~ and At., the preceding tables provide an easy means of checking the 




possibility of E being greater than 100 percent for a given problem. For the case pre- 

K 1 

viously discussed, K.. = 2. 0 and K„ = 0. 5; therefore, - = 4. 0. For At. = 0. 25, the 

1 3 K 3 1 


table indicates T 0 _ = 1. 6. 

2, e, °° 


r (X=0) = 2. 0. 


Since U = 2.0, the tables indicate T„ 

2, max 

Therefore, the possibility exists for an overall effectiveness greater than 100 percent, 
and there will eventually be a decrease in E for an increase in size. 


APPLICATION OF PERFORMANCE CURVES FOR DESIGN 


A complete graphical presentation of performance data is not practical because of 
the large number of independent heat- exchanger variables. The approach used in this 
study is to obtain performance data for selected values of the variables K^, Kg, U, 

and AL which bracket the range of practical interest and to develop interpolation tech- 
niques for intermediate values. While the data presented by the curves are limited, a 
fundamental understanding of three-fluid, crossflow heat exchangers may be obtained 
from them. 

Three sample problems are solved below to demonstrate the application of the 
performance curves for design and to illustrate the physical significance of certain 
trends in the performance data which contribute to an understanding of the performance 
of three-fluid, crossflow heat exchangers. 


Problem 1 

This problem will illustrate the use of the performance curves when no interpola- 
tion is required. It is desired to predict the outlet temperatures of three fluids for a 
heat exchanger operating at the following conditions. 


Fluid 

m, 

lb/hr 

c , 

P „ 

Btu/lb/ F 

V 

h. 

1 

250 

0. 5 

300 

2 

500 

. 5 

100 

3 

250 

1.0 

500 


The surface conductances are 


Uj ? = 50 Btu/hr/ft^/°F and Ug g = 25 Btu/hr/ft‘V°F, 
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while the area x y 
cr o 


• C X4.2 

is 5 ft . 


The resulting nondimensional independent variables are 


m.,c = 

K, = .* P’ 1 - 

1 m 2 c p, 2 


= 0.5 


m„C 9 

= 3 P> 3 = 1>0 

3 m 2 C p,2 


t= . - t_ . > 

At. = =0.5 

1 *3,1 ’Vi 


> (43) 


u„ 


1 9 

U = r— ~ = 2. 0 


u 


2, 3 


u. „x v 

NTU, - l 2 =2.0 
1 m l C p, 1 


Referring to figure 22, the resulting nondimensional dependent variables are 

= 48 percent, 0 ^ = 30 percent, and E = 51 percent. Average outlet temperatures 

t= and t„ may be obtained using equations (44) and (45). 
j.. 6m o • 6m 



1, i 2, i 


j.’ em y x 100 percent 


(44) 


0„ = ( — - m ) x 100 percent (45) 

3 V % i " % i J 

Thus, for 0= = 48 percent, t= = 204° F, and for 0 Q = 30 percent, t„ = 380° F. 
’ 1 1 ’ 1, em 3 c ’ 3, em 

The average exit value of t ^ may be obtained from an energy balance 


m.,c 1 (t 1 . - t= ) + m„c „/t„ . - t„ ) = rii 0 c - t„ 

1 p, 1^ 1, l 1, eny 3 p, 3^ 3, l 3, em/ 2 p, 2^ 2, em 2, l) 


} (46) 


125(300 - 204) + 250(500 - 380) = 250ft - 100) 

i & y 6m j 


thus giving 1 2 em = 268° F. 

It is interesting to note the effect of increasing the size of the heat exchanger to 

2 

x Q y o = 10 ft or NTUj = 4. 0. From figure 22, this change in area results in the fol- 
lowing effectiveness factors: 0^ = 46 percent, 0^ = 44 percent, and E = 66 percent. 
Both E and 0 ^ have increased; however, 0^ has decreased by 2 percent. This 
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1.0 


Fluid (1) .7 


Fluid (3) 


• \ ’ '1 * , 
* ’■ 4 * • 

Added area 

f 

Original 

area 



. . 7 

1.0 


Fluid (2) 

X 




decrease occurs because the increased size of the heat exchanger allows fluid (3) 
(hottest) to heat the center fluid (2) to a temperature greater than that of fluid (1) in 
some parts of the enlarged heat exchanger. 

The dotted region in the sketch indicates the part of the heat exchanger in which 
fluid (2) is hotter than fluid (1). Since fluid (1) is being heated in the dotted region, an 
increase in size actually decreases the effectiveness of fluid (1). Fluid (1) is cooled 
for a smaller value of the coordinate Y. However, as fluid (2) moves through the ex- 
changer, it is heated by both fluids (1) and (3). In the dotted region, heat transferred 
from the hotter fluid (3) increases the temperature of fluid (2) above that of fluid (1). 

The detailed temperature distributions for fluids (1), (2), and (3) are presented 
in figures 49(a), 49(b), and 49(c) respectively. Figure 49(a) indicates that the tempera- 
ture of fluid (1) is above its inlet value of 300° F in almost one-third of the heat 
exchanger. For values of Y less than 0. 45, fluid (1) is always cooled as it flows in 
the X direction. For larger values of Y, fluid (1) is first heated and then cooled as 
it flows through the exchanger. The regions in which it is being heated correspond to 
the regions in which fluid (2), the "coolant" fluid, is actually hotter than fluid (1). 

Figure 49(b) shows that fluid (2) flowing in the Y direction is heated at a very 
high rate along the Y-axis where the hot fluids enter the exchanger. The cooling rate 
decreases for larger values of X. 

As indicated in figure 49(c), fluid (3) flowing in the X direction is cooled most 
rapidly along the X-axis where the coolant fluid (2) enters the exchanger. Since fluid (2) 
is heated as it passes through the exchanger, its ability to cool the outer fluids is de- 
creased. This is illustrated by the isothermal lines for fluid (3) which indicate a de- 
creasing cooling rate for increasing values of Y. 
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In figure 50, isothermal contours for all three fluids are superimposed. This fig- 
ure may be used to determine the part of the heat exchanger in which the temperature of 
the coolant fluid (2) is actually greater than the temperature of fluid (1). This area cor- 
responds to the dotted region. 

Figures 49(a), 49(b), 49(c), and 50 graphically illustrate the cause of the reduc- 
tion of effectiveness for an increase in size which may occur in some three-fluid, 
crossflow heat exchangers. 


This phenomenon is most pronounced for a small At. ^large tg .^, a large capac- 
ity rate for fluid (3) ^large K^, and a small conductance ratio ^large Ug ^ . These 


factors contribute to a high heat-transfer rate from fluid (3) to fluid (2), with the pos- 
sible result that fluid (1) is reheated in some part of the heat exchanger. In some 
cases, for example, K^ = 1.0, Kg = 1.0, and At. = 0.25 (figs. 29 and 31), the effect 

is so pronounced that fluid (1) is actually heated, and 0^ becomes negative as NTU 


increases. The preceding discussion has assumed an original intent to cool fluids (1) 
and (3) with fluid (2). 


Problem 2 


The conditions of problem 2 are chosen to illustrate how the performance curves 
may be used when the independent variables are not equal to those chosen for preparing 
the curves, namely, the combinations resulting from the values listed in the set of in- 
dependent parameters presented in the section entitled "Single-Pass Results. " 


It is desired to determine the temperature effectiveness 9 ^ and 0g for a heat 
exchanger operating under the following conditions: Kj = 0. 40, Kg = 0. 75, U = 1. 3, 
At. = 0.85, and NTU. = 2.0. 

Since these values of the independent variables do not correspond to those for 
which the performance curves have been prepared, some interpolation scheme must be 
employed to determine temperature effectiveness for this case. A straightforward, 
graphical technique has been used. Figure 51 presents the temperature effectiveness 
($ 1 , 63 ) as a function of the conductance ratio U for NTU^ = 2. 0 and At = 0. 85. 

Nine curves are required to cover all possible combinations of K^ and Kg. The data 

were obtained from the performance charts using visual interpolation for At = 0. 85. 

Points were obtained for the three values of U used in the performance curves, 
namely, 0.5, 1.0, and 2.0. Figure 51 was used to determine 9 ^ and 0g as functions 

for Kg for U=1.3 and K^ = 0. 25, 0.50, and 1.0. The results are plotted in parts 

(a), (b), and (c) of figure 52. These curves were used to determine 0^ and 0g as 

functions of Kj for U = 1.3 and Kg = 0.75. Figure 52(d) may be used to determine 

0^ and 0g for Kj = 0.40, Kg = 0.75, U = 1. 3, and At = 0.85. Entering figure 52(d) 
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at = 0. 40 gives 6 ^ = 58 percent and 9 3 = 41 percent. The overall effectiveness 

E could have been obtained in a similar manner; however, it is easier to return to the 
basic definition and calculate it as 


Q a = Kj(l - 



(47) 


or 


Q a = Kjflj + K 3 At.6 3 = 0. 40(0. 58) + 0. 75(0. 85)(0. 41) = 0. 492 
Since + K g > 1 



K i + K a 



K 3 \ _ 1. 28 _ 
+ At.; - 1. is - 


1 . 11 


(48) 


(49) 


Q a 

therefore E = = 44 percent. It is evident from examining figures 51 and 52 that lin- 

' c *oo 

ear interpolation is not adequate; hence, at least three points are needed on each curve. 
An additional point is available on curves which present 0^ and 9 ^ as a function of 

either K- or K„ (fig. 52). Whenever K. approaches zero, 9 . approaches unity. 

3 3 3 

Physically, this means that as the flow rate of a fluid becomes infinitely small, it can ' 
be cooled or heated very easily. It should also be noted that as U becomes very large, 
0 3 will approach zero in figure 51. This trend reflects the physical consequence of 

u 2 3 approaching zero. Insulation of fluid (2) from fluid (3) will not alter the tempera- 
ture. At first it may seem that 0^ should approach zero as U becomes small; how- 


could no longer be 2 . 0 as presumed 


x o y o u i g 

ever, if U approaches zero, NTU 1 = — ? — 

X IT1 ^ C 4 

l p, l 

in the problem. Therefore, it is impossible to consider U approaching zero for a 
finite, constant value of NTU^ . 


When solving a problem where conditions do not correspond to those used in the 
performance curves, a set of figures similar to figures 51 and 52 must be prepared 
for every set of values for NTU^ and At. 
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Problem 3 


Problems 1 and 2 involved predicting output conditions for a given heat exchanger 
operating at specific input conditions. Problem 3 is one which is more frequently en- 
countered by a designer: If the inlet conditions and capacity rates of the two outer 
fluids (1) a.nd (3) are given, determine the size of the exchanger and the mass flow rate 
of the center fluid (2) that are required to produce specified outlet conditions for 
fluids (1) and (3). 

Consider fluids (1) and (3) entering the exchanger at the following conditions. 


riij = 250 lb/hr 
c , = 1. 0 Btu/lb/°F 

p, 1 

t, . = 300° F 


rhg = 250 lb/hr 

c „ = 0. 5 Btu/lb/° F ) 
p, o 

t„ . = 500° F 
3, l / 


(50) 


Coolant fluid (2) is available at T 0 . = 100° F with c „ = 0. 5 Btu/lb/°F. Determine 

i Pj ^ 

the NTUj and m 2 required to cool both fluids to 220 F. These temperature changes 
correspond to the following values of effectiveness. 


9 


1 




300 - 220 
300 - 100 


= 40 percent 


(51) 


and 


9 


3 




500 - 220 
500 - 100 


= 70 percent 


(52) 


K 1 

Possible solutions may be found by studying the design curves for which — - = 2. 0. 

K 3 

This condition is satisfied for two sets of curves: figures 15 to 17 for which = 0. 5 

and Kg = 0. 25 and figures 26 to 28 for which = 1. 0 and Kg = 0. 5. For each of 

these sets of curves, a table may be prepared to investigate possible solutions. First, 
the effectiveness 9 ^ is fixed at 40 percent and the number of transfer units NTU^ is 
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ii mu in ii 


determined for each value of the conductance ratio U. The values of 0 g at this value 
of NTUj are tabulated. A similar procedure is followed holding 0g fixed at 70 per- 
cent and determining 0.. The object is to find a combination of U and NTU ^ for 
which 0 1 and 9 g are as close to the desired values as possible. The following table 
is derived for this problem, for K., = 1.0 and K„ = 0. 5 (figs. 26 to 28). 


u 

0 1 = 40 percent 

0g = 70 percent 

NTUj 

0 3 

NTUj 

0 1 

0. 5 

2. 5 

79 

0.60 

21 

1.0 

2. 5 

75 

1.0 

30 

2. 0 

2. 5 

66 

3. 5 

43 


For U = 2. 0, a value of NTUj between 2. 5 and 3. 5 should be acceptable. Figure 28 
shows that for U = 2. 0 and NTU^ = 3. 0, 9 ^ = 42 percent and 0g = 68 percent, which 
is close enough for design purposes. For this condition the overall effectiveness E is 

riil c j 

83 percent. The required value of m 0 c „ is — — = 250 Btu/hr/°F. Since 

& Pj a 

c 0 = 0. 5 Btu/lb/°F, the required m 9 is 500 lb/hr. 

p, / *2 

It is still necessary to consider the case of = 0. 5 and Kg = 0. 25. A similar 
table is prepared for Kj = 0. 5 and Kg = 0. 25 (figs. 14 to 16). 


U 

0 ^ = 40 percent 

0 g = 70 percent 

NTUj 

0 3 

NTU 1 

9 1 

0. 5 

0. 85 

85 

0. 40 

23 

1.0 

o 

00 

• 

70 

.80 

40 

2.0 

.75 

45 

1.80 

57 


33 




In this case, the values for U = 1. 0 give the desired result exactly for NTU^ = 0. 80. 
From figure 15 the overall effectiveness E is 55 percent. The required value of riig 
is 1000 lb/hr. 

The two sets of conditions which satisfy the objectives of the problem are 


NTUj = 3. 0 
rhg = 500 lb/hr 

> 

U = 2.0 

E = 83 percent / 


(53) 


and 


NTUj = 0. 80 

rhg = 1000 lb/hr 

> 

U = 1.0 

E = 55 percent / 


(54) 


The choice facing the designer is between a large physical size with low flow rate and 
high effectiveness or small size, larger flow rate, and lower effectiveness. The ulti- 
mate choice must be based on factors such as construction, cost, space available, vol- 
ume of coolant fluid available, and other design factors. 

Not all problems which are approached in this manner will have an adequate solu- 
tion. Consider a case which has the same conditions as the previous problem except 

K 1 

that ttt - has the value 0. 50 instead of 2. 0. That is, assume that the hotter fluid has 
K 3 

the higher capacity rate. Performance curves must now be considered for which 
Kj = 0. 25 and Kg = 0. 50 (figs. 8 to 10), and Kj = 0. 50 and Kg = 1. 0 (figs. 20 to 22). 

A chart is prepared as before for K^ = 0. 25 and Kg = 0. 50. 
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u 

9 1 = 40 percent 

0g = 70 percent 

NTUj 

*3 

NT^ 

0 1 

0. 5 

0.85 

48 

2.0 

53 

1.0 

.70 

26 

4.0 

67 

a 2. 0 

— 

— 

— 

— 


cL 

0 j is always above 9 ^ for U = 2. 0. 


The trends indicate that there will be no satisfactory solution. For K^ = 0. 5 and 

Kg = 1. 0, 0g never even reaches 70 percent. In this case, the performance curves 

indicate that the objectives of the problem are impossible under the imposed restric- 
tions. 


CONCLUDING REMARKS 


Performance characteristics have been determined for a wide range of operating 
parameters for single-pass, three-fluid, crossflow heat exchangers. The performance 
of two-pass heat exchangers for both cocurrent and countercurrent flow has been stud- 
ied for selected operating conditions. The results have been presented in terms of the 
temperature effectiveness of the two outer fluids as functions of heat- exchanger size 
for sets of fixed operating conditions. 

Selected values have been chosen to bracket the range of practical interest be- 
cause of the infinite possibilities for combinations of operating conditions. Interpola- 
tion techniques have been used to obtain performance data for intermediate values. 
Sample problems are included to illustrate the use of the performance curves and the 
interpolation techniques. 

An expression for overall effectiveness has been derived which compares the heat 
transferred by a particular exchanger with that which is transferred by one of infinite 
size. Isolated cases corresponding to poor design are cited for which the overall effec- 
tiveness may be greater than unity. This indicates the importance of using the temper- 
ature effectiveness of the two outer fluids as the primary design variables and the 
overall effectiveness as an auxiliary parameter. 

While data are necessarily limited to fixed sets of operating conditions, a funda- 
mental understanding of three-fluid, crossflow heat exchangers may be obtained from 
the performance curves. 
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A computer program has been developed for the study of both single- and 
multiple-pass heat exchangers. Output options are available for detailed studies of 
temperature distributions within a particular exchanger and for generation of perform- 
ance data for a large number of heat exchangers. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, November 17, 1967 
905-89-00-00-72 



APPENDIX A 


NUMERICAL PROCEDURE 


Basic Logic 

To determine the temperature distributions of each fluid in a three -fluid, cross - 
flow heat exchanger, the partial differential equations which must be solved simulta- 
neously are 




(Al) 


0T„ 

- W =C { T 1- T 2) + D ( T 3~ T 2, ) < A2 > 




(A3) 


The region of solution of these nondimensionalized equations is the portion of the 
X-Y plane bounded by X = 0, X = 1, Y = 0, and Y = 1. The boundary conditions are 
T 1 = T 1 (Y) and T 3 = T 3 (y ) at x = 0 and t 2 = T 2 (X) at Y = 0. 



0 t 2 = T 2 (X) 1 X 
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The basic logic for solution will be outlined before the details of the integration 
scheme are discussed. 

1. The numerical integration of equations (Al) and (A3) could be initiated in the 
X direction if Tg were known along the Y-axis. 

2. To obtain Tg(0, Y), equation (A2) is integrated numerically in the Y direc- 
tion. 

3. Using the initial values of T and Tg and the values of Tg at X = 0 calcu- 
lated in step 2, T and T g can be calculated at X = AX using equations (Al) and 
(A3). 

4. At X = AX, the same situation exists as before: T and T„ are known and 

1 o 

Tg is to be calculated from equation (A2). 

5. The above procedure is repeated at each increment AX until the solution is 
obtained over the entire region. The only difference in the initial integration step and 
all the other steps is that T and Tg are no longer constant along an X = constant 

line. 


Integration Scheme 

The numerical technique used was devised as a first-order, predictor -corrector 
integration scheme. The solution of equation (A2) for Tg(Y) at X = X N will be used 

to illustrate the procedure. Assume that T , T„, and T_ are known at X = 

Y = Y n and the value of Tg is desired at Y = Y^ + AY. It will be recalled, from the 

outline of the basic logic, that T and T„ are known along the line X = X from 

13 N 

( dT 2\ 

Y = 0 to Y = 1. Equation (A2) is used to evaluate -J at Y= Y^. This derivative 

/* t 2 \ 

will be denoted by ^ prediction of the value for Tg at Y = Y^ + AY is 

Y N 

calculated from 




+ 



(A4) 
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This value of T 2 P is then used in equation (A2) along with the known values of 

/9T 2 \ 

Tj and at Y = Y N + AY to predict l -gyi at Y = Y N + AY. This derivative will 


' 3T, 


be denoted by l J 


Y n+ AY 


. A corrected value of T 2 = T 2 is calculated from 



(A5) 


This procedure can be shown to be equivalent to using a second order Taylor 
series expansion of the function T 2 (Y) at the point Y N with the required first deriv- 

3Tj 9T 3 

atives and approximated by the slopes between Y N and Y^ + AY. 

A similar procedure is used to solve equations (Al) and (A3) in the X direction. 
In the solution of these equations, T 2 is assumed to be constant over the interval AX 

between X^ and X^ + AX^ at the value T^X^. 


Accuracy Check 

The accuracy of the computation may be checked at any x coordinate during the 
integration by comparing the energy gained (lost) by fluid (2) with that which is lost 
(gained) by fluids (1) and (3). 

The calculation of the energy balance proceeds in the following manner: 
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At any station x N> tj, and t g are averaged from y = 0 to y = y Q while 1 2 is aver- 
aged from x = 0 to x = x N> Conservation of energy requires 



The coordinates — may be replaced by X N> Dividing equation (A6) by /m 2 c A 
x o VP’/ 

the resulting equation is 


vv*! 2 ) 




+ 


A 3 C p, 3\ 

A 2 C p, 2 ] 




(A7) 


1T1 j C w o 

Since . 1 p ’ = K. and — — = k„, equation (A2) may be written 

m 2 C p, 2 1 m 2 C p, 2 3 


m„c 


K 


lQi, i " t l, m( X N)] + i " t 3, m( X N)] X l|\ m( X N) " fc 2, i] (A8) 


Dividing by ^ i " l 2 i gives 


K l[ T l, i - T l, m( X N)] + K 3 [ T 3, i - T 3, m ( X N)] = X n[ T 2, m( X N )J < A9 > 


Since the boundary conditions are . = 1 and the accuracy of the 

overall computation may be checked at any station X N by comparing the quantities 

K l[i - T l, rn^)] + K 3 aE" T 3, mN 31111 X n[ T 2, m( X N )} ™ s <=»“Pari S on is 
used to determine the appropriate step size for the different calculations. 
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An examination of the basic equations 



would disclose the direct influence of A, B, C, and D on the calculation. 

For large A, B, C, or D, the temperature gradients will be large, and a 
smaller step size will be required to maintain acceptable accuracy. Since any of these 
constants is a nondimensional representation of the size of the heat exchanger, larger 
exchangers will require more calculational steps. 

It was arbitrarily decided that an acceptable limit for accuracy would be that the 
two overall energy -balance terms would not differ from each other by more than 2 per- 
cent of their average value. 

An automatic step -size control was used in the computer program because of the 
large number of cases which were needed to generate the performance curves. For 
any calculation, the largest value of the set (A,B,C,D) is denoted as M. The follow- 
ing criteria for step size were established. 


M 

Ax, Ay 

M >20 

0. 002 

7 < M <20 

.005 

4 < M <7 

.01 

2 < M <4 

.02 

1 < M <2 

.05 

M < 1 

. 10 
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In all cases, this set of criteria was sufficient to insure agreement of the energy 
balance within 2 percent, and the agreement was considerably better in the majority of 
cases. 
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APPENDIX B 


COMPUTER PROGRAM 


A computer program has been developed to solve the basic differential equations 
using the numerical procedure described in appendix A. The program is capable of 
handling calculations for both single- and two-pass heat exchangers. 

Since the program is written in FORTRAN, it can be run on the IBM 7094 or 
Univac 1107/1108 computers. The following pages contain a listing of the complete 
program. 
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PROGRAM MAIN 


DIMENSION Tl(lOOl), T2CG01), T3 Cl 001 > , XAVE (1 4 ) , XI I (1 001 > , Y 1 1 (100 
XI) , BCDX (12) • BCDYC12), T2PCX00X). T1 T ClO) , TIP (501 > , T3P<501) 

1 , THPPP C 1 001 ) 

COMMON /AB/ A. B, C, 0, DTI 
REAL A 1 , A3 
C 
C 

A COO = D 

C 

READ (5(16) (XAVE(AJ), AJ=1,14 ) 

DO 40 JO=l ,14 

40 XAVE(JO) = XAVE(JO) -.001 

I r ( XAVEC1) .LT. 0.00001 .AND. XAVEC2) .LT. 0. 00001 ) XAVE ( 1 > = 1.0 
I ZERO = 0 
CALL RESET 
100 CONTINUE 

RE... (5,995) IPASS, I TYPE , IHIX, IDENT, IPLOT, IPRINT, I ALL 
3722 CONTINUE 


DTI, A, U 
(1 .0/U> 


kEAU(5,IS) A1 , A3, 

C = A * A1 
B = A * (A 1/7.3) * 

D = (Al/U) * A 
WRITE <6, 1000) 

IP < I TYPE .EQ. DWRITE (6, 1C01) 

IF ( I TYPE .EQ. 2JWRITE (6, 1010) IPASS 
IF < I TYPE .EQ. 3)WRITE<6, 1011) IPASS 


101 


IF <1 TYPE 

.EQ. 

. 1 ) GO TO 101 



IF (IKIX 

.EQ. 

0) WRITE (6, 1012) 



ifcimix 

.EQ. 

1 ) WR I TE (6,1013) 



IPCIOENT 

.EQ. 

, 0)WRITE(6, 1014) 



ifcident 

.EQ. 

. 1 ) WR ITE(6,l0l5) 



CONTINUE 





A3C0 = A 





IFCA3CD 

.LT. 

B) ASC0=B 



ifcaocd 

.LT. 

C) A3CD=C 



IFCABC0 

.LT. 

D) A3C0=D 



NX = 500 





IP (ASCD 

.LE. 

20.0 .AN'O. AOCD .GE 

. 7.0) 

NX = 200 

IP ( A3CD 

.LT. 

7.0 .AND. ABCO .GE. 

4.0) 

NX - ICO 

IP (ASCD 

.LT. 

4.0 .AND. ABCO .GT. 

2.0) 

NX = 50 

IP (ABCO 

-LE. 

2.0 .AND. ABCO .GE. 

1.0) 

NX = 20 

IP (ABCO 

.LT. 

1.0) NX = m 



X 

Z 

It 

>- 

Z 





NN = NY 

♦ 1 




N?1 = NX 

+ 1 




OELX = 1 

.0 / 

FLOAT (NX) 



dely = 1 

.0 / 

FLOAT (NY) 




WR 1 TE (6 , 1006) A , A1 
WRITE (6, 1009) B , C, 
WRITEC6, 1016)DELX, 
306 CONTINUE 
Til □ 1.0 
722 = 0.0 


A3 , 
D 

OELY 


U, DTI 
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T33 = 1.0 / OTI 

T2GUES = 0 . 5+ ( (K1+K3/DTI)/ (K1+K3) ) 

UTx I 7 £ C 6 , 1 00 7) T3 3 

IFUTYPE .£Q. 3) VH ITE(6,2S) T2GUES 
IF { I TYPE .EG . 3>T22 = T2GUES 
IF 1 1 PLOT .cO. 0) GO TO 304 
XII (1) = 0.0 
YI J (11 = 0.0 
CO 41 IK = 1 t HP 1 
XIIUK+l) = XII (IK) + DELX 
YIMIK + I) = YII(IK) ♦ OELY 
41 CONTINUE 
304 CONTINUE 
IAVE = 1 

IF ( I PASS .EQ. 1) IAVE = 0 
WHITE (6, 17) 

INV = 0 

DO 42 166=1 ,10 
177 = 1 
203 CONTINUE 

IF < I NV .EQ. 1 ) GO TO 200 
DO 1 1=1 ,NN 
T1 (I) = Til 
T2P ( I ) = T22 
l T3 ( I ) = T33 
200 X = 0.0 
INV = 0 
KK = 1 

T2SUN = 0.0 
IXO = HOO ( I 77 , 2 ) 

I F ( I PASS .EG. 1 ) GO TO 3 
WH I T c (6,29) 166 
I F ( I NO ,N2. 0)V/3ITE(6,21> 

I F ( I NO .EC. 0) ViH I T£ (6 , 22) 

3 DO 4 J = 1 , HP 1 

CALL AB1 (Tl, T2, T3 , DELY, NY , T2P(J) ) 
T2P(J) = T2(NN) 

T2SUH = T2SUM T2 <NN) 

KCOU = KCCXJ ♦ 1 

IF (KCGU .NE. 1 PR 1 NT) GO TO 202 

KCOO = 0 

WH!TE(6,11)X 

W3IT£(6,12) 

V3 I T£ (6 , 10) (Tl <K) ( K=1,NN> 

W^ITE(6, 13) 

WH ITE(6, 10) (T2(K) ,K=1 ,NN ) 

I Tl (6, 14) 

VHITEC6.10) <T3 (K) ,K=1 ,NN ) 

202 CONTINUE 

I F ( XAVE(MO ,LT. O.DGO TO 371 
IF(XAVE(AK.) .LT. 0. IE-03) GO TO 18 
IF (XAVE (KK) .GT. X) GO TO 18 
IF < I ALL . EQ • 0) GO TO 310 
371 CONTINUE 

WRI 7E(6, 11)X 



WRXTc (6.221 

WT?I7£(6, 10) (T1 (K),K=1,NN) 
W.?!TECo,13> 

W.?I7£<6.10) (T2 (K) ,K = 1 ,NN ) 

WHITE <6, 14) 

WHITE(6,10) (T3 (K> ,K = 1 ,NN ) 

31 G KK - KK ♦ 1 

T2S = 72SUM / FLOAT (J) 

TtSUH = 0.0 

T3SUM = 0.0 

DO 19 IK=1.KN 

T 1 SUM = 71 SUM ♦ T1CIK) 

19 T3SUX = T3SUH ♦ T3(IK) 

II SUM = T1SUH / FLOAT (NN) 
T3SUM = T3SUH / FLOAT (/AN) 
IFCPLOT .SsZ. l>CO TO 302 


CALL 

cu:k:»:l 

(-1 . 

o 

o 

1.0, 

a 

o 

5.0, 

1H1 , DCOX, BCDY , NN , Y I I ,T1) 

CALL 

CUZKKL 

( o* 

0.0, 

1 .0, 

0.0, 

5.0, 

IK 2 , DCOX , OCOY ,NN , Y I I ,T2) 

CALL 

GUI XML 

( o. 

o 

o 

1.0, 

0.0, 

5.0, 

1H3 , OCDX , BCDY ,NN , YI I , T3) 


302 continue: 

W3 I TE (6 , 1 G2 1 ) 

WHITE(6,20> T1SUH, T2S, T3SUM 

CON = K1 * (1.0 - T 1 SUN) ♦ K3 * (1.0/DTI - T3SUK) 

XX22 = 0.0 

Ir ( I TYPe .EQ. 2>XX22=T22 
CO.M1 = X * (T2S-XX22) 

IF( IKO .NZ. 0) COM 1 = (T2S ~ T22J * X 

IF( I WO .EQ. 0 ) CON = XI * (Til - T1SUH) ♦ K3*(T33 - 

1 T3SUH) 

ACC = ADS (CGH-COH1 ) / (COM ♦CONI ) * 1GO.O 

IF ( I TYPE .LE. 2)WHITE(6, 1004) COM, CW1, ACC 
18 CONTINUE 

X = X ♦ CZLX 
TIP ( J) = T 1 (NN) 

T3P ( J ) = T3 (NN) 

CO 2 1=1 ,NN 

DEL = DELX * A * (T2(I) - T1(I) > 

T1P1= T1 (I ) ♦ DEL 

DELI = DZLX « A * (T2(I) - TIPI) 

Tl ( I ) = (DEL + 0EL1HC.5 ♦ T1(I) 

DEL = CZLX * B * (T2 (I ) - T3(I) ) 

Tl?l= T3 (I) ♦ CEL 

DELI = DELX * 6 * (T2<I)-T1P1 ) 

2 T3 ( I ) = 73(1) ♦ (DZL+0CL1) * 0.5 
4 CONTINUE 

CALL CUZKKL (-1 ,0.0, l .0,0. 0,5.0, 1H1 ,CC0X ,ECDY ,N.N, Yl I , TIP ) 
CALL CUZKKL ( 0 , 0 . G , I . 0 , 0 . 0 , 5 . 0 , 1 H? , DCOX , BCDY , NN , Y I I , T2P ) 
CALL 0 U I KML ( 0 ,0.0 , 1 .0,0.0, 5.0, !H3 , DCOX , BCDY ,NN, YI I , T3P ) 

X = 1 .G 
311 CONTINUE 


TriEl 

= 1 .0 

- 

Tl SUM 

THE 2 

= DTI 

* 

T2S 

TKE3 

it 

a 

- 

DTI 4-T3SUH 

C\ = 

K 1 ♦ 

K 3 


CAVE 

= 0.5 

V 

(CON + COH1 ) 
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C'i A* = (IU * K3/uT«) * (iiO / GU 

IF(GA .LT. 1.0) Q-rl AX = K1 ♦ K3/DTI 
E = CAVE / CMAX 
V3I7E<6, 1005) 

3000 CONTINUE 

IF(I>*.0 .NE. 0) T2S3S = T2S 
IF ( I AVE .EQ. 0)GO TO 100 
177 = 177 ♦ 1 
IF < I 77 .GT. 2) GO TO <4 

1FC1XXX -NE- TO <6 

711 = TiSUH 
T22 = 0.0 
T33 = T3SUM 

IFUTYPc .EQ. 2) T22=T2S 
CO TO 203 

45 CONTINUE 
T22 * T2S 
Til = 1.0 

T33 = 1.0 / DTI 
IF ( 2 HI X .EQ. 1) INV = 1 
IFUDENT .EC . 0) GO TO 52 
INV = 1 
DO 9 155=1 ,NN 
T 1 (155) = Til 
T3 (155) = T33 
9 THPPPCI55) = T2P(I55) 

DO 91 155=1 , NN 
1551 = NN - 155 ♦ 1 
91 T2PCI55) = TMPPP (15 51) 

GO TO 42 

46 CONTINUE 

IF (IDE NT .EQ. 0) GO TO 47 
DO- 4S 155 = 1 • NN 
49 TMPPP (15 5) = T 1 C I 55) 

DO 5 155=1 , NN 
1551 = NN - 155 ♦ 1 

5 T 1 (155) = TMPPP (1551) 

DO 6 155 = 1 i NN 

6 TMPPP (155) = T3C255) 

DO 7 15 5 = 1, NN 

1551 = NN - 155 ♦ 1 

7 T3 ( I 55) = TMPPP (1551) 

DO & 155=1 , NN 

8 T2P(I55) = 0.0 
Til = T 1 SUM 
T22 = T2S 

T33 = T 3 SUH 

1F< 1 TYPE .EQ. 2) GO TO 56 
GO TO 200 

56 DO 57 155 = 1 , NN 
1551 = NN - 155 ♦ 1 

57 T2P(I55) = T2U551) 

GO TO 200 

52 CO 53 NKO=l,NN 
Tl (KKO) = T 1 1 
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53 T3(KX0) = T33 

IF< I K] X .EG. 1) INV = 

CO TO 42 

47 CONTINUE 
IF Cl TYPE .EQ. 2) CO TO 58 
DO 43 155=1 f NN 

48 T2P < I 55) = 0.0 
CO TO 200 

58 DO 59 KKO=l,NN 

59 T2PCNK0) = T2(KKO> 

CO TO 200 

44 CONTINUE 

IF ( I TYPE .L£. 2) CO TO 103 

Tl T ( I 66) = T2S 

IF ( 166 .EQ. l)CO TO 45 

IF ( (ASS ( Tl T (1 66- 1 ) - Tl T (I 66) ) /Tl T (I 66- 1 ) > .LT. 0.01 ) CO TO 103 

Til = 1.0 

T22 = T2S 

T33 =1.0/ DTI 

IF (I DENT .EQ. 0)CO TO 52 

INV = 1 

DO 102 155=1, NN 
Tl (155) = 1.0 
T3U55) = T33 

102 TNPPPCI55) = T2P(I55) 

DO 92 155 = 1 , NN 

1551 = NN - 155 ♦ 1 
92 T2P (155) = TNFPPU551) 

42 CONTINUE 
WR IT£(6,23) 

103 CONTINUE 
W3ITE(6, 1022) 

THElO = 1.0 - T 1 SUM 
THE30 =1.0- 0TI * T3SUM 
020 □ T2SSS 

IF ( I TYPE .EQ. 2) G20 = T2S 

C 1 30 = M * (1.0-T1SUM) * K3 * Cl. 0/DTI - T3SUH) 

EO = (0130 ♦ 020) * 0.5 / CMAX 

ACC = A5S (Q20- Q 1 30) / (Q130+Q20) * 100.0 

QAVEO = (02CH-C13O) * 0.5 

WRITE(6, 1004)0130, 020, ACC 

I F (O*^ .LT. 1 .OJW'RITE <6, 1003) 

IF (OX .CE. 1.0) W3 I TE (6 , 1002) 

AA2 = A ♦ A 

WRI 7£ (6, 1024) CAVEO, QHAX , AA2 , EO, THElO, THE30 
CALL TIME Cl TIME) 

I ZERO = HIKE - I ZERO 
U^ITc (6, 1029) I ZERO, I TIME 
I ZERO - I Tl HE 

1029 FORMAT (1H0, 10X, 19HTIHE FOR THIS CASE ,I5,13H MICROSECONDS, // 

1 1 1 X , 1 1 HTOT AL TIME ,I8,13H H I CROSECONDS ) 

W3I7E (6, 1005) 

CALL OKPCUF 
CO TO 100 

10 FORHAT ( 6 ( 5X , E 1 4 . 7) ) 



11 FORMAT<1HO,///,5X,4HX = ,F5.3 ) 

12 FCF<NAT<1H0,///,2X,16HT1 (I) , I = 1,N ) 

13 FORMAT <1H0,///,2X,16HT2(I) , I = 1,N ) 

14 FORMAT ( 1 HO , /// , 2X , 1 6HT3 ( I ) , X = 1 ,N ) 

15 FCT.; !AT (214, 611 ) 

1G FORMAT (7r 1G.2) 

17 FORMAT (1H0, //, 55X,14H>** OUTPUT *** ) 

20 FORMAT (1H0 ,40X , 26HTHE AVERAOc value OF Tl = ,F8.3, / 41 X t 

1 2SHTHE AVERAGE VALUE CF T2 = ,F8.3, / 41X, 

2 26HTKE AVERAGE VALUE CF T3 = ,F8.3 ) 

21 FOR . .7<1N0,57X, 8KCOMAIN I ) 

22 FORMAT (1H0.57X, 9HDCMAIN II ) 

23 FORMAT (1HQ,2QX,45HTHE SYSTEN DID NOT CONVERAGE IN 10 ITERATIONS) 

28 FCRMAT(1H0,46X,24HT2 CUESS FOR DOMAIN I = ,F5.2) 

2D FORMAT ( 1 KO . 5GX , 10HI TER AT ICN ,12 ) 

55 FORMAT (12A6/12A6) 

SS5 FORMAT ( 7211 ) 

1CC0 FORMAT ( 1 H 1 , 53X,13K^^ INPUT > 

1C01 FORMAT < 1H0 , 59X , 1 1H5 1 NGLE PASS ) 

1002 rCRMAT(lrtQ,20X,3S:-iAl A3 IS GREATER THAN OR EQUAL TO 1.0 > 

1GG3 FORMAT <1HQi20X,2*JKK1 A3 IS LESS THAN 1.0 ) 

1004 FORMAT { 1 H . 35X, 44HFGR ACCURACY CHECK ON ENERGY BALANCE COMPARE// 

1 2SX, 6KC13 = , F8. 3, 13H WITH Q2 = , F8 . 3 , 5X , 3HOR ,F7.4 i 

2 1 SH PER CENT ACCURACY ) 

1C05 FORMAT < 1 HO , 128 ( 1 H>) ) 

1007 FORMAT (1H0, 40X, 23H0GUNDAR Y CONDITIONS / 43X • 

1 SHT1 = 1.00,5X,8HT2 = 0.0.5X.5HT3 = ,F4.2 ) 

1003 FCRMAT(lHu,59X, 12HINPUT VALUES / 30X,4KA = , F4 . 2 » 5X , 5HK 1 = , 

1 r 4 . 2 , 5 X , 5HK3 = ,F4.2,5X,4HU = , F4 . 2 , 5X , 1 1 HDELTA TI = ,F4.2 ) 

1009 FORMAT (1H0 , 40X, 3SHTHC RESULTING VALUES OF B, C, AND 0 ARE / 

1 4 3 X , 4H3 = ,F5.2, 5X, 4HC = ,F5.2, 5X, 4HD = ,F5.2 ) 

1010 FORMAT <1 HO, , 51X, 1 5KMULT I PLE PASS (.12. 10H) PARALLEL FLOW) 

1011 FORMAT (1H0, , 51X, 1 5HMUL T I PLE PASS (,I2, 17H) COUNTER FLOW ) 

1012 FOR MAT ( 1 H , 6RX , 5HMIXEO) 

1013 FORMAT <1H , 6) X i 7HUMMIXE0 ) 

1014 FORMAT (1H , 57X, 1 4 M I N VER TED ORDER ) 

1015 FC.ViAKlH , 5 7 X , 15HI0ENTICAL ORDER ) 

1016 FORM AT ( 1 r.G , 50X , lOhDELTA X = ,F5.3 / 5 1 . X » 1 CjHDELTA Y = ,F5.3 ) 

1017 FOR MAT (lH0,40X t 45HDATA IS PRINTED OUT AT CACH CALCULATED POINT. ) 
1013 FCRMAT(lHu,40X,46KDATA IS PRINTED CUT AT THE FOLLOWING VALUES OF X 

1) 

1019 FORMAT! 5 1 X , 2HX ( , 1 2 , 4H> = , Fj.3 ) ^ 

1020 FCr ; :jAT(lH0,40X,7H2AVE - ,F8.3 , 10X.7HCMAX = ,F8.3 // 50X, 

1 4HA = , F 8 . 3 //30X,4hc = ,F8.3 , 1 OX , 9HTHET A 1 = ,F8.3 ,10X 
2, CHTHETA3 = ,F6.3 ) 

1021 FOR MAT (1H3,34X,33HTHE AVERAGE EXIT TEMPERATURES ARE ) 

1022 FORMAT(lHO,41X,32r{THE OVERALL OUTPUT VARIA3LES ARE) 

1023 FCRrUT (1H0.35X ,44H"C.\ ACCURACY CHECK ON ENERGY BALANCE COMPARE // 

1 2SX, 7H013O = ,F0.3,14H WITH QZO = ,F8.3 ,5X, 3HOR ,F7.4 , 

2 18.-; PER CENT ACCURACY ) 


1024 FORMAT ( 1 HO , 3 3X , 3i-iG A VEO = ,F8.3, 10X , 7HQMAX = ,F8.3 // 48X, 

17KN7U1 = ,F6. 3//24X.5HEO = , F8 . 3 , 1 OX , 1 CHTHETA 1 O = ,F8. 3.10X, 

2 1GHTHETA30 = , F8.3 ) 

END 



SUBROUTINE ABi <Ti , T2 , T3 , OELY , NX , T2& > 

COMMON /AB/ A, 6, C, D, DTI 
DIMENSION T1 (1001 > » T2<1001>, T3<1001> 

T2<1> - T20 
OO 1 N=1 , NX 

DEL = DELY * <C*T1 <N> ♦ 0*T3 <N) - (CO) * 72 (N> 

TEKPT2 = T2<N) ♦ DEL 

DELI = DELY * <C*T1 CN*1) ♦ D*T3<N+1> - <OD> * 

1 T2(N+l> = T2(N) ♦ (DEL+0EL1) * 0.5 
RETURN 
ENO 


) 

TENPT2 ) 
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Figure 1. - Schematic representation of a single -pass, three -fluid 
heat exchanger in crossflow. 
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Figure 2. - Schematic representation of the heat -transfer surface 
between fluids (1) and (2). 
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Figure 3. - Schematic representation of a countercurrent two-pass, three-fluid, crossflow heat exchanger. 
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Figure 7. - Effectiveness factors for a single -i 

K 1 = 0.25, K = 0. 25,' 


Overall effectiveness, E, percent 
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U = 2. 0. 
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Figure 9. - Effectiveness factors f< 
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Figure 13. - Effectiveness factors for a single -pass, three -fluid, crossflow heat exchanger; 

K 1 = 0. 25, K- = 1.0, and U = 2.0. 
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Figure 16. - Effectiveness factors for a single -pass, three -fluid, crossflow heat exchanger; 

K 1 = 0. 50, K„ = 0. 25, and U = 2. 0. 
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Figure 20. - Effectiveness factors for a single -pass, three -fluid, crossflow heat exchanger; 

K 1 = 0. 50, K„ = 1. 0, and U = 0. 50. 
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Figure 30. - Effectiveness factors for a single -pass, three -fluid, crossflow heat exchanger; 

= 1.0, K„ = 1.0, and U = 1.0. 
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Figure 32. - Effectiveness factors for a two -pass, countercurrent, three -fluid, crossflow heat exchanger 
with mixed flow in elbows; K 1 = 0. 50, K„ = 0. 50, and U = 0. 50. 
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Figure 33. - Effectiveness factors for a two -pass, countercurrent, three -fluid, crossflow heat exchanger 
with mixed flow in elbows; K 1 = 0. 50, K„ = 0. 50, and U = 1. 0. 
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Figure 34. - Effectiveness factors for a two-pass, countercurrent, three-fluid, crossflow heat exchanger 
with mixed flow in elbows; K., = 0. 50, = 0. 50, and U = 2. 0. 
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Figure 36. - Effectiveness factors for a two -pass, countercurrent, three -fluid, crossflow heat exchanger 
with identical flow order in elbows; K 1 = 0. 50, K„ = 0. 50, and U = 1. 0. 
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Figure 37. - Effectiveness factors for a two -pass, cc 
with inverted flow order in elbo\ 
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Figure 38. - Comparison of effectiveness factors for single -pass and two-pass heat exchangers; 

K 1 = 0. 50, K„ = 0. 50, U = 1. 0, and At. = 0. 25. 



Figure 39. - Convergence of Tg ^ for pass 1 in a two-pass, countercurrent heat exchanger 

with inverted flow in the elbows. 
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Figure 40. - Temperature distributions in two -fluid, counterflow exchangers 
with finite and infinite heat -transfer area. 
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(a) Fluid (1). 

Figure 49. - Isothermal contours. 
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(c) Fluid (3). 
Figure 49.- Concluded. 



i 

t 

* 

-S-67-6776 


lu id 

m, 

Ib/hr 

Cn, 

Btu/lb/ °F 

°F 


!1) 

250 

0.5 

300 

Ui # 2 = 50 Btu/hr/ft 2 / 0 F 

:2) 

500 

.5 

100 

U2,3 = 25 Btu/hr/ft 2 / 0 F 

3) 

250 

1.0 

500 

x„v„ = 10 ft 2 
























National Aeronautics and Space Administration 
Washington, D. C. 20546 

official business FIRST CLASS MAIL 


POSTAGE AND FEES PAl4 
NATIONAL AERONAUTICS A 
SPACE ADMINISTRATION? 


08U 001 58 51 3DS 68120 00903 p 
AIR FORCE WEAPONS LABORATORY/ AFWL / [ 
MRTLAND air FORCE BASE, NEW MEXICO 871H 


ATT MISS MADELINE F e CANDY A , CHIEF TtCHNI,. 
LIBRARY / W L I L / i' 


If 


POSTMASTER: 



If Undeliverable ( Section 13^ 
Postal Manual ) Do Not ReftT 



"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the ividest practicable attd appropriate dissemination 
of information concerning its activities and the results thereof.” 


— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


I 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


